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Abstract: A food web on n living things 71, 72,---,%n, i.e., a biological n-system 
can be mathematically characterized by action flow Cc. of order n with surplus flows 
of growth rates #; of population on vertices v;, vector flow (z;,x;), end-operators 
rifij, cj fj, on edge (v;,vj), where fi;, fj; are 2 variable functions for integers 1 < 


1,9 <n holding with a system of conservation equations 
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which is a system of n differential equations. Certainly, 0 € R” is one of its equilibrium 


points. But the system 
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of equations may be solvable or not. However, even if it is non-solvable, it charac- 
terizes biological systems also if it can be classified into solvable subsystems. The 
main purpose of this paper is to characterize the biological behavior of such sys- 
tems with global stability by a combinatorial approach, i.e., establish the relationship 
between solvable subsystems of a biological n-system with Eulerian subgraphs of la- 
beling bi-digraph of (om characterize n-system with linear growth rate and the global 
stability on subgraphs, and interpret also the biological behavior of G’-solutions of 
non-solvable equations, which opened a way for characterizing biological system with 
species more than 3, i.e., mathematical combinatorics. As we know, nearly all papers 
discussed biological system with species less or equal to 3 in the past decades. 
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§1. Introduction 


There is a well-known biological law for living things in the natural world, i.e., the 
survival of the fittest in the natural selection because of the limited resources of 
foods. Thus, foods naturally result in connection with living things, i.e., food chain, 
a linear network starting from producer organisms and ending at apex predator 
species or decomposer species. And biologically, a food web is a natural interconnec- 
tion of food chains, a resultant by a simple ruler ([28]), and generally a graphical 
representation of what-eats-what in the ecological community such as those shown in 
Fig.1 for 4 food chains: grass—ladybug—frog—snake—eagle, grass—ladybug—frog 


—eegret, grass—rabbit— snake—eagle and grass—rabbit—eagle. 


Fig.1 


Actually, a food web is an interaction system in physics ([15]-[16], [25]) which 
can be mathematically characterized by the strength of what action on what. For a 
biological 2-system, let x, y be the two species with the action strength F’(x — y), 
F(y > x) of x to y and y to x on their growth rate, respectively ({21]). Then, such 
a system can be quantitatively characterized by differential equations 
g£=F(y— 2) 
y= F(z y) 
on the populations of species x and y. 
Usually, we denote 2 competing things by a directed edge (u,v) labeling with 


vector flow (x,y) and end-operators F’, F” respectively on its center and both ends, 


where F’, F’ are action operators with F(x — 0) = F’(0 > y) =0ify=O0orzr=0 
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and the growth rates £, y of populations on vertices, such as those shown in Fig.2. 
Particularly, F = «f, F’ = yf’ in the Kolmogorov model, where f, f’ are 2 variable 
functions, and f = — by, f’ = w+ cx in the Lotka-Volterra model ([2], [20]). 


G) F (x,y) F! ® 


Fig.2 


Then, a food web is nothing else but a topological digraph ren a 2-tuple (v(G), E(G)) 
with E(G) C V(G)xV(G) anda labeling L: G — RU SonG with L: V(G) > 
Rand E(G) > S, where R and S are predetermined sets ({19]). Particularly, if 
R= {x,y}, the growth rates of populations and S = {(F, (x,y), F’)}, a 3-tuple with 
action operator F’ on the initial, F’ on the end and vector (x,y) on the middle of 
edge (u,v), we get the biological 2-system shown in Fig.2. 

However, the law of conservation of matter concludes that matter is neither 
created nor destroyed in chemical reactions. In other words, the mass of any one 
element at the beginning of a reaction will equal to that of element at the end, i.e., 

: : 2 : aL 
the in and out-action must be conservative with the surplus on each vertex of G . 

D 
Thus, a food web is an action flow ([18]) further, ie., a topological digraph @ 
labeled with surplus flows of growth rates 4; of population on vertices v;, vector flow 


(x;,7;), initial and end operators Fi, Fy; on edge (v;,v,;) for integers 1 < i,7 <n, 


where n > 2 holding with a system of conservation equations 
t= SS F(t, > xi) — SS Fi(a;— 2), 1<i<n 
UREN (vi) vwEN* (v4) 


and particularly, 


B=a)| > Seene=~ >, -faleyn) |, 1 Stee". Ga) 
vEEN (1%) vwEN* (v4) 
in the Kolmogorov model. For example, a biological 4-system shown in Fig.3 is a 


system of 4 ordinary differential equations 


1 = 21 (b52 — a1 — Qy2)£1 — A12%2 — A22%3 + bs124) 


Lo = Ly (bi121 + (bie + beg — a31)%2 — A32%3 + 6124) 


(1.2) 


boi % + b31%2 + (b22 + b32 — G41) %3 — A42%4) 


b41X3 — 452%, — Ag2X2 + (ba2 — a51 — ag1) £4) 


t3 = © 


——_~ ~~ 


La = U4 


3 


where, 
fi(@1, 22) = An%1 + Q12%2, fi (@1,%2) = but + Oi2%2, 
fo(x1, 3) = G10 1 42203, fp(1, r3) = by121 + bo2%3, 
f3(@2, £3) = A312 + 3273, £3 (2, £3) = b31%2 + bg2x3, 
fa(xs, 4) = 44103 T G42X4, fi (3, 4) = bai 13 + bazXa, 
fs(@4, 01) = G51%4 + G52%1, f$(@4, 21) = b51%4 + b5271, 
fe(@4,%2) = Ag1%4 + Ag2%2, fE(@4, Lo) = be1a4 + be2Xo. 


Definition 1.1 Let GC be a labeling topological digraph. A subgraph H of G is said 
to be a labeling subgraph of Ga if its vertices and edges are labeled by L|,,, denoted 


by ie < eG and furthermore, if Ae = Cilhes such a labeling subgraph is said 


to be an induced subgraph of (els denoted by (V(F)) 


For example, the 2 labeling graphs Gi, G in Fig.4 are all labeling subgraphs 
but only G, is an induced subgraph of the graph shown in Fig.3. 


Fig.4 


Clearly, a labeling subgraph of Ge is also consisting of food chains but it maybe 


not a food web if it is not an action flow again. Even it is, the sizes of species are 


i 
not the same as they in G” because the conservative laws are completely changed. 


i 
For example, the system of conservation equations for the labeling subgraph oe is 


Ly, = 2 ((b51 = 91) X4 — 92%3 + bs14) 
£3 = 3 (ba 21 + (bag — G41) %3 — G42%4) (1.3) 
La = XL4 (ba x3 — A52%1 + (baz = 51) 4) 

a very different system from that of (1.2). 


The following terminologies are useful for characterizing food webs. 


Definition 1.2 Let G be a digraph with G a digraph reversing direction on ev- 
ery edge in G.a bi-digraph of G is defined by GiuLe and a labeling bi-digraph 
feue) of a labeling digraph Ge ws labeling graph on Giuie with a la- 
beling L : V(@UG) = L(v(@)), L: E(GUG) + L(z(GUG)) by 
L: (u,v) > 10, (2,9), 9J", (v,u) > {xf,(z,y), OF PL: (u,v) > {2f, (ey), yft 
for V(u,v) € E(G), such as those shown in Fig.5. 


Cig) 
rf (vy) yf’ ae : _ 
Glen wp n& oO 
x 0 
(x,y) 


Fig.5 


Definition 1.3. A circuit in a digraph G is a nontrivial closed trail with different 

edges in G and an Eulerian circuit in digraph G is a circuit of red containing every 
—=> 

edge of G. 


A digraph G is Eulerian if it contains an Eulerian circuit. 


Clearly, a bi-digraph of a digraph is an Eulerian graph. The main purpose 
of this paper is to characterize the biological behavior of biological n-systems with 
global stability by a combinatorial approach, i.e., establish the relationship between 
solvable subsystems of a biological n-system with that of labeling Eulerian subgraphs 
of labeling bi-digraph (Gus) of @: characterize conditions of an n-system 
with linear growth rate become distinct and global stability, and interpret also the 


biological behavior of G’-solutions of non-solvable equations, which opened a way 


for characterizing biological system with species more than 3, i.e., mathematical 
combinatorics, or differential equations over graphs. 

For terminologies and notations not mentioned here, we follow references [1] for 
mechanics, [25] for interaction particles, [2] and [20] for biological mathematics, [3] 
for differential equations with stability, [6]-[7] for topological graphs, digraphs and 


combinatorial geometry, [7] and [26] for Smarandache multispaces. 


§2. Geometry Over Equilibrium Points 


2.1 Equilibrium Sets 


We consider the generalized Kolmogorov model on biological n-system ([2], [20]), 


i.e., the system (1.1) of differential equations 


satisfying conditions following: 


(1) fy, fj € C’ for integers 1 < i,j <n; 
(2) For any integer 7, 1 <i <n, there is (x?, 78,---, 2°) € R” hold with 


an 


VEEN (v4) vpENt (vj) 


but 


S- Of 
( k? i) vpENt (vj) 


aft, 
rs 


UREN (vj) 


For any integer 7, 1 <i <n, define 


F,= S° Fiei(Lk, Bi) — S- Fil Ze 2p- 


vREN (%) yENt (vi) 


Then it concludes that 


(1) F; € C! for integers 1 <i <n; 
(2) Fi] (29 28,29) = 0 but 


Ofne 
= SS r : 


x 
(x9 ,29,---,29) UREN (vi) 


OF; 
OX; 


Se 


x 
(x9 ,29,---,29) vpEN* (vj) 4 


#0. 


(x?.29,--2h) 


Applying the implicity function theorem, each equation 
Ft tis gio: _ “4 ie) = 0 


is solvable, i.e., there is a solution manifold Sp, in R” for any integer 1 <i<n, 


and in this case furthermore, there is a unique solution on the Cauchy problem 


of system (1.1) prescribed with an initial condition (x1(to), V2(to),-++,¢n(to)) = 
(zt, i, one) t): 

An equilibrium set of system (1.1) are all points (r?,79,---, 2°) € R” holding 
with 


tt ( De fea (2h, a) = » f(z, “) =0 


UREN (v1) vpEN* (v1) 


a XS foleta)- ¥ fale) =0 (2.1) 


trl dD  fenl@e%n)—- dF  fuilen, 7) | = 0 
UREN (Un) yENt (vn) 


Clearly, only those solutions x? > 0,1 <i <n of system (2.1) have the biolog- 


ical meaning, and (0,0,---,0) € IR” is an obvious equilibrium point. We classify all 


equilibrium points of system (2.1) into 3 categories following: 


(C1) Only (0,0,---,0) € R” hold with system (2.1), i.e., the system 


» fia (eh, 21) = ye f(x}, 2?) 


UREN (v1) vtEN (v1) 
» fho(ah, £5) = > fala, x?) 
UREN (v2) upENTt (v2) (2.2) 


ee case = me FD ster) 


UpEN- (vn) yENt (vn) 


is non-solvable in R”. 


(C2) Only (0,---,0, A1,0,--+,0, Ko,0,---,0, Ks,0,---,0) € R” hold system 


(2.1) with numbers Ky, Ko,---,K; > 0 on columns 71, i2,---,i, respectively, i.e., 


for any integer j ¢ {i1,72,---,i,}, the system 


3 fat.) a iS Fit(2},, 22) 


UREN (vi, ) wEN* (vi, ) 
2.3 
FE Phalehah) = S falah a) ee 
UREN Vis yENT Vis 
2D fag (Uh, 24) = Ss fix, 2?) 
UREN — (v;) vpENT (v5) 


is non-solvable in R”. 
(C3) There are (Ay, Ko,---, Ky) € R” hold system (2.1) with K; > 0 for 
integers 1 <2 <n. 


2.2 Geometry Over Equations 


Usually, one applies differential equations to characterize the reality of things by 
their solutions. But can this notion describes the all behavior of things? Certainly 


not ({8]-[19]), particularly in biology follows by the discussion following. 


For an integer n > 1, let u: R” — R” be differentiable mapping. Its n- 
dimensional graph I'|u] is defined by the ordered pairs 


Plu] = eae 7 -, En) »U (Lig on “y))) (eas a gaa) 


in R"*!. Clearly, the assumption on f;;, fij,1 < 4,7 <n concludes that the solution 


manifold Sr, is nothing else but a graph [[F;] in R”. 


Geometrically, the system 


Fy (M89, 4594) = 0 
Fo(21, £2,°-+, Ln) =0 (2.4) 
Geen = 

is solvable or not dependent on 0 Sr, 0 or not, and conversely, if . Sr, #9 or 
i= i=1 

not, we can or can not choose point (21,%2,+-++,%p) in a Sr, a solution of (2.4). 


i=1 
We therefore get a simple but meaningful conclusion following. 
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Theorem 2.1 A system of equations 


under previous assumption is non-solvable or not if and only if () Sp =0 or #0. 
i=1 


If the intersection () Spr, £0, it is said to be a A-solution of equations (2.4). 


i=1 

Usually, one characterizes a system S of things 7), 7>,---, 7}, by equations (2.4) 
with their solutions to hold on the dynamical behavior of these things. Js zt always 
right? The answer is negative at least in the non-solvable case of equations (2.4), 
and even if they are solvable, it can be used only to characterize those of coherent 
behaviors of things in S, not the individual such as those of discussions on multiverse 
of particles in [15] and [16]. Then, what is its basis in philosophy? It results deeply 
in an assumption on things, i.e., the behavior of things discussed is always consistent, 

e., the system (2.4) is solvable. If it holds, the behavior of these things then can 
be completely characterized by the intersection Q Spr,, Le., the solution of system 
(2.4). However, this is a wrong aigieteandine on things because all things are 
in contradiction in the nature even for human ourselves, and further on different 
species. This fact also concludes that characterizing things by solvable system (2.4) 
of equations is only part, not the global, and with no conclusion if it is non-solvable 
in classical meaning. 

Philosophically, things 7), 7>,---, 7}, consist of a group, or a union set U T;, and 
if T; is characterized by the ith equation in (2.4), they are geometrically equivalent 
to the union U Sr,, e., a Smarandache multispace, not the intersection a Sr. 

For oa le: if things T), T>, 73,74 and Ty, T;, 73, Tj are respectively Tanita 


ized by systems of equations following 


x+y = 1 Lr=y 
+ = —| =") 
(GHSNy a 8 (pss) 4 774 
e=ap = S1 r=1 
Pays = 


then it is clear that (LESJ’) is non-solvable because x + y = —1 is contradictious 
to x+y = 1, and so that for equations x — y = —1 and x — y = 1, ie., there are 
no solutions Zo, yo hold with this system. But (LES?) is solvable with 2 = 1 and 
y =1. Can we conclude that things Ti, T5,T3,T, are x = 1,y = 1 and T,, To, T3, Ty 
are nothing? Certainly not because (x,y) = (1,1) is the intersection of straight line 


behavior of things T7, 73,73, 7; and there are no intersection of T;, 7, 73, 7 in plane 


R?. However, they are indeed exist in R? such as those shown in Fig.6. 


Fig.6 


Denoted by the point set 
Lape = {(@,y)|ax + by = c,ab # 0} 


in R?. Then, we are easily know the straight line behaviors of T,, 7»,73,74 and 
Ti, 13,T3,T; are nothings else but the unions LZ) -19U Lii2U Li01U Loi, and 
Lyag U Ly 3-1 U | be U Tit. cg respectively. 


Definition 2.2 A \-solution, also called G-solution of system (2.4) is a labeling 
graph G” defined by 


ViG) ={Sm, bax ak 
E(G) = {(Sm, Sp,) if Se,(\Se, #0 for integers 1 < i,j <n} with a labeling 
L: Sr Sr, (Sr, Sr,) > Sp, ()Sz,. 


Example 2.3 The V-solutions of (LES!) and (LES?) are respectively labeling 
graphs C/ and Kf shown in Fig.7 following. 
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Fig.7 


Theorem 2.4 A system (2.4) of equations is V-solvable if F; € C' and Fi] (29 28,29) 
OF; 


Oi | (49 29,09) 


= 0 but # 0 for any integeri, 1<i<n. 


Proof Applying the implicity function theorem, the proof is completed by 


definition. 

iz 
Theorem 2.5 A system (1.1) of differential equations on food web G’ is uniquely 
V-solvable if fiz, fi; € C’ for integers 1 < i,j <n and (a1(0), #2(0),--+,4n(0)) = 
(x9, 29,---, 2°) ER”. 


Proof Applying the existence and uniqueness theorem on the Cauchy problem 


of differential equations, 


Ly = Xj y Fit (Sieve) S Fula, £1) 


VRE N- (w i vpEN* (vj) 


with (21(0), 72(0),--+,%n(0)) = (xo, 7$,---,2°) € R”, it is uniquely solvable for any 


integer 1 <i <n. Consequently, the system (1.1) is uniquely V-solvable in R” by 


definition. 


2.3 Equilibrium Sets of Linear Equations 


Certainly, the Lotka-Volterra model on biological 2-system is a system of linear 
growth rates. Generally, if all fj;, fj; are linear for integers 1 < 7,7 <n, then it 
is a generalization of Lotka- Volterra model on biological n-system. We can further 


characterize the equilibrium sets of linear system (2.4) by linear algebra. 
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Definition 2.6 For any positive integers i,j, i # j, the linear equations 
QD + O3F2b Cindy = Oj, 
G27 - Ggabe Ff ** ja tn = 0; 
are called parallel if there exists a constant c such that 
C= GG = Gea] Ge HS SS Gp) Gin 0; 10; 
The following criterion is known in [8]. 


Theorem 2.7([{8]) For any integers i,j, i # j, the linear equation system 


Gi Fy Optet aig te = bj, 


Aj1L1 + AzoL2 +--+ AjnTn = b; 
is non-solvable if and only if they are parallel. 


By Theorem 2.7, we divide all linear equations L;, 1 < i < n in (2.4) into 
parallel families 


C1, 62,°°+, Gs 


by the property that all equations in a family @; are parallel and there are no 
other equations parallel to equations in @; for integers 1 < i < s. Denoted by 
\G;| = ni, 1 <i < s. Then we can characterize equilibrium sets of linear system 
(2.1) by Theorem 2.6 in [8] following. 


Theorem 2.8(|8]) The equilibrium sets of system (2.1) with linear growth rates 


fii, fig, VS t,9 <n can be classified into 3 classes following: 


(LC1)_ there is only point (0,0,---,0) € R” holding with linear system (2.1), 
1.e., its V-solution 


Ce kK 


M1 .NQ,";Ns 


with ny +n+2+4+---+n, =n, where n; = |G;| and G; is the parallel family for 
integers 1<i<s,s>2. 

(LC2) there is only point (0,---,0,c1,0,--+,0,¢2,0,--+,0,Cn_1,0,---,0) € R” 
holding system (2.1) with numbers cy, C2,-++,Cn—1 > 0 respectively on columns 14, io, 
+++ ,dn_ for 1 <I <n, i.e., its V-solution 


i E 
G = ie tee 
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with nw t+n+24+---+m = 1, where n; = |G| and G; is the parallel family for 
integers 1<i<t, s>2. 
(LC3) there is an unique point (C1, C2,+++,€n) € R” holding linear system (2.1) 


with constant c; > 0 for integers 1 <i <n. 


§3. Biology Over Equations 


Classically, a solvable system (1.1) of differential equations characterizes dynamical 
behaviors of a food web in area. However, the solvable systems are individual but 
non-solvable systems are universal. Then what about biology over those of non- 
solvable systems (1.1)? Are there no biological significance? The answer is negative. 

Firstly, let us think about a food web how to run. Certainly, a food chain only 
follows a direct, linear pathway of one animal at a time, and different thing T’ has 


his own food chain for living, even for the same kind of things. 


Fig.8 


For example, the eagle can lives respectively on the rabbit, on the snake or on 
the both via its food chains snake—eagle or rabbit—eagle with interactions in Fig.1, 
i.e., although the eagle preys on the snake and the rabbit but it is also dependent 
on the 2 populations such as those shown in Fig.8, and its living web should be 
consisted of circuits eagle—snake—eagle, eagle—rabbit—eagle or eagle—snake and 
rabbit—eagle, Eulerian subgraphs. 

Generally, a predator P preys on a living thing 7’, i.e., P action on T and there 
are also T’ reacts on P at the same time, which implies the interaction between 
living things, the in and out action exist together. We therefore know a biological 


fact following. 


Fact 3.1 A living thing must live in an Eulerian subgraph of bi-digraph of a food 


web Ga’. 
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The following result characterizes action flows on Eulerian subgraphs with that 


of solvable subsystems of equations (1.1). 
Theorem 3.2 Let G. be a food web with solvable or non-solvable conservative 
equations (1.1) on initial value (x1(0), v2(0), +++, @n(0)) = (x9, x8,--+, 2°) and A» 


L L 
(G U G) , a food web containing species T with solvable conservation equations 


Lig = Lig ys Fig 2k, Lio) = S- Figt (Lio, £1) », Lsts || (3.2) 


UREN (vig) WEN (vig ) 


ii 
in H where L(0;.)- = Sig. Then H is an Eulerian digraph and H® is an action 


flow. 


Proof If He is a food web, by Fact 3.1 H must be an Eulerian digraph. 

Now let x, = f(%1, £2,°++,2%n),1 <i < |H| be the solution of (2.5). Notice 
that in solution x;,, x; can be any chosen constant c, particularly, 7; = 0 if 7 ¢ 
4105. UE S| i (32) es 


Lig =F (Ql 30; wag. 0 0 te, 032s Og Oe Oe a | | 


is also the solution of (2.5) with fix(vi, 7%) = 0, fi;(an, vi) = 0, which implies that 


L, F : ; 
H is an action flow with conservation laws at each vertex. 


Let Ay, Ho, ee, H, be subgraphs of digraph G with H, H; ie NA; = 

or # @ for integers 1 <4 j < s. If Gz U Ai, they are called a subgraph multi- 
i=1 

decomposition of Ga Furthermore, if each H; is Eulearian, such a decomposition is 


called an Eulerian multi-decomposition, denoted by C= B A. For example, an 


i=1 
Eulerian multi-decomposition of the graph on the left is shown on the right in Fig.9. 
e V1 U1 U1 
U3 
— v3 @ v4 DB Us 
U4 U5 
V2 (0D) U2 
v2 


Fig.9 
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Particularly, if E (1;) (VE (7;) = () for integers 1 <i #7 < s, such a decompo- 
sition on G is called an Fulearian decomposition of G. 


The next result characterizes food webs by Eulerian multi-decomposition. 


‘s 
Theorem 3.3 I/f there are Eulerian subgraphs H,; , L<i<s with solvable conser- 


Dee, hee 
vative equations, i.e., food webs such that G = DB H, with 
i=1 


Liv ae Yu e€ VIG) 


l 
ifue f) V(Hin) with L(v) = Zi in H,, and 


ip=1 


(sR uv) Con SA w), V(u,v) € E(G’) 


[=1 


if (uv)e 1) E(H), then G is also a food web, i.e., an action flow. 
jo=1 
Proof Clearly, em is a labeling graph holding with conservative law on each 


aL. , 2 
vertex u€ V(G ), i.e., an action flow. 


84. Global Stability and Extinction 


In biology, the generation is the necessary condition for the continuation of species 
in a food web constraint on the interaction, i.e., the stability in dynamics with small 
perturbations on initial values. Usually, the dynamical behavior is characterized by 
differential equations, which maybe solvable or not and can not immediately apply 
to the stability of food web a for n > 3 by Theorems 3.2 and 3.3. Generalizing 
the classical stability enables one to define the stability of food web following. 


Definition 4.1 A food web @” with initial value Ee where L(v) = £y, Lo(v) = £9 
forve V(@’) is globally stable or asymptotically stable for any initial value Ge: 
where Ly(v) = y2 forv € vie ) and a number €, > 0, there is always a number 
dy > 0 such that if \y° — x°| < 6, exists for allt > 0, then 


y(t) —ae(t)]<e,, WoeV(G), 
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or furthermore, 
SL 
lim |y_(t) —2,(t)|=0, WeEV(G ). 


Certainly, we need new criterions on the classic for discussing the stability of 


species in biology. 


Theorem 4.2 A food web em with initial value a is globally stable or asymptot- 


ically stable if and only if there is an Eulerian multi-decomposition 
LE . L 
(¢Ua)'-O7 
i=1 
with solvable stable or asymptotically stable conservative equations on labeling Eule- 
L 
rian subgraphs H, for integers 1<i<s. 
9 iD 
Proof The necessary is obvious because if G with initial value G ° is globally 


iB 
stable or asymptotically stable, then (G U G) is Eulerian itself by Fact 3.1. 


Now if there is an Eulerian multi-decomposition 
fz =. =Si 
i=1 


re 
on labeling bi-digraph (G U G) with stable or asymptotically stable conservative 


equations on labeling Eulerian subgraphs H i.e., for any number ¢, > 0, there is 


1) 


a number 6, > 0 such that if |y? — x°| < 6, exists for all t > 0, then 
£ 
lyo(t) — x,(t)|< ey, Woe V(A;), 
or furthermore, 


lim |y(t) —2,(t)| =0, Wo €V(H,) 


t0 


i 
for integers 1 <2 < s, let A, be the multiple of vertex v € ges ) appeared in 
L 
subgraphs H, ,l <i <.s, we then know that 


lyo(t) — 2,(t)| < ef 


L 
for v € V(E; ) if jy? — x°| < 6% for integers 1 <i < Ay. 


7 


Define 
6,=min{di,1<i<A,} and e,=maxf{et,1<i<A,}. 
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We therefore know that 

lyo(t) — to(t)| < ev, 
i.e., the species on vertex v is stable if the conservative equations of H. are stable 
for integers 1 <2 < A, and Ge is globally stable. 


Now if furthermore, x, is asymptotically stable, i.e., 


lim |yu(t) > EAt)| = 0 


t0 


ie 
in food web H, ,L<i< Aj, it is clear that 


lim lyv(t) 7% Ly(t)| = 0 


£ i 
in G also, i.e., ag globally asymptotically stable. This completes the proof. 


Corollary 4.3 A food web Ga with initial value ae is globally stable or asymp- 


totically stable if there is an Eulerian decomposition 
(@U2) -@F: 
i=1 
with solvable stable or asymptotically stable conservative equations on labeling Eule- 
rian subgraphs ries for integers 1 <i<s. 
Clearly, the bi-digraph eal) G has an Eulerian decomposition, called parallel 


GUG = OB ((u.v) Ue.) 


—=> 
(uv)EE(G) 


decomposition 


We get the next conclusion. 


Corollary 4.4 A food web Ga with initial value Ge is globally stable or asymp- 
totically stable if it is parallel stable or asymptotically stable, t.e., ((u,v) Uv, u))” 
is an action flow for V(u,v) € E(@’). 

For an equilibrium point G@ et (2.1), we can also linearize F'(v,u), F’(v, u) at 
(0, Yo) for V(v,u) € E(@’) and know the stable behavior of Gin neighborhood 


of ec by applying the following well-known result. 


Theorem 4.5([3]) Jf an n-dimensional system X= F(X) has an equilibrium point 


Xo that is hyperbolic, t.e., all of the eigenvalues of DF, have nonzero real parts, 
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then the nonlinear flow is conjugate to the flow of the linearized system in a neigh- 
borhood of Xo. 


The next result on the stability of food webs is an immediate application of 
Theorem 4.5. 


Theorem 4.6 A food web @ with initial value Ge is globally asymptotically stable 


if there is an Eulerian multi-decomposition 
B 7 L 
k=1 


with solvable conservative equations such that Red; < 0 for characteristic roots Xj; 

of Ay in the linearization AyXy = On,xn, Of conservative equations at an equilib- 
L is re L 

rium point H,. in H, for integers 1 <i < hy andv € V(A,), where V(A,) = 


{V1, V2,°° "5 Vis 


Qi, Ag *** Ap, 
Ay =| 43, @. +++ Op, 
Ay Aho “°° Oh, 
a constant matric and Xp = (Xy,,Lv.5°°*,Ly,,)" for integers 1<k <1. 


Proof Applying the theory of linear ordinary differential equations, we are 


easily know the species 
hy 
B= Ga ae", 
i=1 


where, c; is a constant, 3;(¢) is an h,-dimensional vector consisting of polynomials 


in t determined as follows 


ti tyit + ti 
= toi = tait + too 
y(t) oe i Bol ) ad | ene : 
thy tait + th,2 
Toes a TT eo tik, 
t bat t ky -2 
ice (t) = ie es at ; Pak 
1 9 
Get + tant? Choks 


7 to(k1+1) > _ tat + toe 
By as) ’ Br, a(t) ~e ’ 
Uhy (ki +1) tnit + th,2 
SE eee eset Oh es Nien tl tis ater Nt Len ee Bl acta ; 
H(hy—kst1) pks—1 ) '(ny—kst2) ¢ks—2 
(ks—1)! t T (ks—2)! t fh 
fa(hy—ks+1) pks—1 4 !2(hy—ks+2) phs—2 ton 


Oe (t) = (ks—1)! (ks—2)! 


thy (hy—kst+) pks—l 1 thu(hy—kst2) pks—2 |. 
a a a t ae S| a t + + thin, 


with each t,; a real number for 1 < 7,7 < h, such that det([tij],, .,,) #9, 


At; tf LS ski; 


Ao, tf ky tl <i<ko; 
Q= 


Age 4)- het e hee e Le tS is 


If Red; < 0 for integers 1 <2 < hy, it is clear that 


lim x,(t) =0 


t—co 


E 

for vertex v € E (F,), ie., each linearized conservative equation A,X, = 0n,xh, is 
L 

stable for 1 << k < s. Applying Theorems 4.2 and 4.5, we therefore know that @ 


is globally asymptotically stable. 


Comparatively, we also get the next conclusion on the unstable of a species by 


Theorem 4.2 following. 


Corollary 4.7 A species T is unstable in a food web rel with initial value a af 


and only if the subgraph containing T in all Eulerian multi-decompositions 
oo GS ue 
i=1 


of ie is unstable. 


A unstable behavior of species T will causes the redistribution of flows and 


makes for a stable situation on the food web Ce If it established, the food web 
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works in order again. Otherwise, a few species will evolve finally to extinction, i.e., 
ceases to exist in that area. If all species in a food web Gr vanished on that area, 
there must be a series of species 71, 72,---, 2; Successively died out on the time, the 
stability of the web is repeatedly broken, established and broken, and finally, all 
species become extinct. In this case there must be vertices v1, V2,---,Us € v(@’) 


and a series of action flows 


=L =L =>L 


G —-G ~14(@ -u)-wo-- 3G — {v1, V2,--+, Us} 


such that there are no flows in @ {U1, U2,°°°,Ust, Le., C= {U1,V2,°°+,Us} & Ki, 
where | = |G —s. 

Notice that if species x is extinct, there must be tim x(t) = 0. Let f(t) be 
a differentiable function on populations of a species x. If f(t) = O(t-*), a > 1, 
i.e., there are constants A > 0 such that |f(t)| < At~° holds with t € (0,+00), 
then f is said to be a-declined and x a species extinct in rate a. Furthermore, if 
f(t) = O(e—) for B > 0, because 


2 
t Bt+ Ft re ae ees eee 


we are easily know that there is a constant A > 0 such that |f(t)| < At” for any 
integer n > 1. In this case, f is said to be oo-declined and x a species extinct in 
rate oo. 


The results following characterize the extinct behavior of species in a food web. 


Theorem 4.8 Let Ga be a food web hold with labeling L : vu; — x; on vertices v;, 
L: (uj, v3) > (Fi, (a, 23), FG} on edges for integers 1<i,7 <n, VC V(@’). If 


VQ=So( SO F's 2)- SS Bw’) 


vEeV \v/EN-(v) v'EN*t(v) 


is a-declined, then all species X = > x, in V is extinct in at least rates a—1 and 
vEV 
particularly, if V = {v}, the species x, is extinct in at least rates a —1 ont. 


Proof Notice that the conservative equation at vertex vu € V is 


C= Ss" F",,(@yt > Ly) — S- Fyy (Ly > Ly) 


v'EN~ (v) v'ENt(v) 
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and 


ip OE ste 


veV 
= S- S- FY,(v' > v) — Ss" Fyy(v > 0’) 
vEeV \v/EN-(v) v'EN*t(v) 


Now, if V(t) is a-declined, there must be constant A > 0 such that 


ci ‘ poe ' A 
vEV v'EN~(v) v'EN+(v) 
Consequently, 
A +00 4 5 
ae t@ (a =, Dee 1 
0 0 


Therefore, all species X in V is extinct in at least rates a — 1 on t, and partic- 
ularly, it holds with the case of V = {v}. 


Theorem 4.9 Let GC be a food web hold with labeling L : uv; — x; on vertices v;, 
Lz (uj, 0;) > {Fiy, (a4, 23), Fi} on edges for integers 1<i,j <n, andV Cc v(G") 
a cut set with components C1, C2,---,C, in C Vi. heres > 2. IT 
A()= SY) Fle’ )- SY) Fw’) 
v'EN~(v) v'EN+(v) 
18 Qy-declined for Vu € V with a = min Qy, then 


1 CG turns to | food webs Ce. Ge . Gk finally; 
ag l 


(2) the species Xy = D> Ly, particularly, x, is extinct in at least rates a — 1 


vEV 
ont for Vu EV. 


Proof Applying Theorem 4.8, all species X in V is extinct in at least rates a—1 
on t, and finally, extinction if t — oo. In this case, there are only left components 
C1, C2,-+-,C, and each of them is a food web because if x, = 0, there must be 
F(v > u) = 0 and F’(u > v) = 0 for Vu € V andu € v(@") \ V. Therefore, the 
conservative laws 


i= S- FY (2u > Vy) — Ls Fie (Sy By) 


vEN-(u) vEN*t(u) 


Pall 


. ob 
in G turns to 


Li a Ss" Bay _ Geis) _ ys Bgl te ars Gi 
veN-(u) 1 V(Ci) veN*(u)(\ V(Ci) 


ra 
i.e., it holds also with vertex u in @. for integers 1 <i <1, the assertion (1). 
For (2), by the proof of Theorem 4.8 there is a number A > 0 such that 


+00 +00 
ae ar A 
— f eG) Sat < Xv(t) =S~ 2,(t) < [@ zat 
0 vEV 0 


by definition, where n(G’) is the connectivity of ce Whence, Xy(t) = O(t-**?), 
and particularly, v,(t) = O(t~**") for v EV. 


Finally, there are indeed the case of extinction of species in rate oo. For example, 


the proof of Theorem 4.6 implies the case of extinction in rate oo on t following. 


Theorem 4.10 Let GC. be food web with an Eulerian multi-decomposition 
L Ss 
(Ua) -@r 
k=1 


z 
and all conservative equations on H, are solvable for integers 1 <k <l. For 
if es i L 
a vertex vu € ViG ) including repeatedly in H, Te Aig dt if ReA; < 0 for 


a4? uy? 


characteristic roots A; of A, in the linearization 


Ay Xe = On tie, 
L L 
of conservative equation at an equilibrium point H,, ve V(H;,) for integers 
1 <2 < hg, then the species xy is simultaneously extinct in rate co on time t 
i 
and asymptotically stable, where V(A,) = {v1, v2,-++, Un, }, Ax and Xz are as the 


same in Theorem 4.6 for integers 1<k <l. 


Proof By the proof of Theorem 4.6, we know that x,(t) is asymptotically stable 
with 


hr 
ee x ci, (t)er* 
i=l 
Define 3 = min A;. If A; < 0 for integers 1 <i < hx, then x,(t) is clearly an 
SUSE 


oo-declined function and species x, is extinct in rate oo on time t. 
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§5. Algorithm 


Let Y = {G, G,---,Gm} be solvable Eulerian multi-decompositions of bi-digraph 
(G U Gy of a food web G’ with conservation equations (1.1) solvable or not, where 
G@, and G,, are respectively a parallel decomposition, rem itself of G. Theorems 
4.2 and 4.6 conclude the following algorithm on the global stability of reds 


Algorithm 5.1 The stability of a food web Ga can be determined by programming 
following: 

STEP 1. Input X; = 6 andi = 1,2,---,m; 

STEP 2. Determine Eulerian circuits in X; is globally stable or not; 

STEP 3. If X; is globally stable, go to STEP 6; Otherwise, go to STEP 4; 

STEP 4. Replace X; by Xj;41, return to STEP 2; 

STEP 5. If Xj41 is globally stable, go to STEP 6; Otherwise, go to STEP 4 if 
i<m, or go to STEP 7 ift =m; 

STEP 6. G@: is globally stable, the algorithm ts terminated; 

STEP 7. a: is globally non-stable, the algorithm is terminated. 


This algorithm certainly enables one to determine the stability of a food web 
L 
G regardless of whether its conservation equations solvable or not, and get stability 


of food webs with more species than 3 by conclusions on 2 or 3 species. 


L 
Example 5.2 Determine the stability of a biological 5-system G” shown in Fig.10, 


Fig.10 
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where, fi; and fj;, 1< i,j <5 are defined by 


for(@2, 21) = 1— 22 — 4X1, fox (42, 01) = i A2Xo, 
fi(a3,%1) =1—2%3-Aia1, — fg,(@3,71) = 1 A3V3, 
fa(t4,41) =1—24-—A1m1, fy (44,41) =1 Naa, 
fal@s, = 1 See Mis. Feiles, ei) = AsZ5 


with conservative equations 


Ly = (4 = 401 = AnL2 = A3X3 = AL = A5i5) 


£2 = —£2 (1 — re — A121) 
£3 = —23 (1 — x3 — A121) (5.1) 
tq = —¥4 (1 — &4 — A121) 
f5 = —£5 (1 — 25 — A421) 


Let (a9, 2$, 7, 29, 2) be an equilibrium point of (5.1). Calculation shows the 


linearization of (5.1) is 


= Ax} = Agxt re = A3x9 x3 = Annan = Asxoxs 


Lo = Ayxhry (-1 Ons d 29) v2 


0 


L3 = Axx, 


(-1 1) 3 
cA = Ayr ry (-1 255 Do ea 
(-1 ) 5 


222 + Ay 29 


Ls = Aizex1 


where A = 4— 829 — Agx$ — A3x$ — Agr} — Asx?. 


As usual, we can hold on the stability of system (5.2) of linear equations and 
then, the stability of (5.1) by Theorem 4.6 on equilibrium points with tedious cal- 
culation. However, we apply Algorithm 5.1 for the objective. 

Notice that bi-digraph ee U GE of @” in Fig.11 has a parallel decomposition 


such as those shown in Fig.12, 


Fig.12 


and the conservation equations on these parallel edges are respectively 


Ly = x1(1 — M1 A2X2) 


Lo = x2(1 — Xa A121) 


twill =e Ages 


t3 = = 23(1 == 23: N21 


ies L— 2, — Agu, 


( 
4 = ta(1 — La — AG Pi 
( 


face: 1— a2, — Ast5 


(5.6) 
v5 = x5(1 om OF ios A121 


) 
) 
) 
(5.5) 
) 
) 
) 
We have known the stability of equations (5.3)-(5.6) by their linearizations 
following ([{20]): 


(1) the equilibrium point (21, x;) = (0,0) is unstable for equations (5.3)-(5.6), 
wheres = 2,3,4;5: 

(2) the equilibrium point (a1, 27;) = (1,0) is stable if Ay > 1 for i = 2,3,4,5; 

(3) the equilibrium point (21, x7;) = (0,1) is stable if A; > 1 for equations (5.3)- 
(5:6), where ¢—=2,3)4.5° 
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A; -1 »4-1 
AA; — VA - 1 
and X; > 1 for equations (5.3)-(5.6), where i = 2,3,4,5. 


(4) the equilibrium point ( ) is asymptotically stable if A; > 1 


Therefore, we know this biological 5-system is unstable on the equilibrium point 
(0,0,0,0,0) but stable on the equilibrium points (0,1,1,1,1) and (1,0,0,0,0), and 


asymptotically stable on the equilibrium point 


MA-VAAH-V AMAH-L AAH- TV AYA 1 


of system (5.1) if A > 1 and A; > 1 by Theorem 4.2, where X\ = Ag = A3 = Ag = As. 


§6. Conclusion 


Today, we have many mathematical theories but we are still helpless on opening the 
mystery of the nature as Einstein’s complaint, i.e., as far as the laws of mathematics 
refer to reality, they are not certain; and as far as they are certain, they do not refer 
to reality because the multiple nature, or contradiction is universal on things, par- 
ticularly, living things different from rigid bodies. Usually, we establish differential 
equations for characterizing things T and holds on their behavior by solutions, which 
is only hold on those of coherent behaviors of things, not the individual. Thus, we 
encounter non-solvable cases in biology, and even if it is solvable, finding the exact 
solution is nearly impossible in most cases. In fact, the solvable equation is indi- 
vidual but the non-solvable is universal for knowing the nature. This fact implies 
that we should also research those of non-solvable equations for revealing reality of 
things in mathematics, which finally brings about the mathematics over topologi- 
cal graphs, i.e., action flows, or mathematical combinatorics, and only which is the 


practicable way for understanding things, particularly, living things in the world. 
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